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ABSTRACT 

A two-dimensional honeycomb lattice harbors a pair of inequivalent valleys 
in the k-space electronic structure, in the vicinities of the vertices of a hexag- 
onal Brillouin zone, K±. It is particularly appealing to exploit this emergent 
degree of freedom of charge carriers, in what is termed " valley tronics" , ^"^ if 
charge carrier imbalance between the valleys can be achieved. The physics 
of valley polarization will make possible electronic devices such as valley fil- 
ter and valley valve, and optoelectronic Hall devices, ^'^ all very promising for 
next-generation electronic and optoelectronic applications. The key challenge 
lies with achieving valley imbalance, of which a convincing demonstration in 
a two-dimensional honeycomb structure remains evasive, while there are only 
a handful of examples for other materials. We show here, using first princi- 
ples calculations, that monolayer M0S2, a novel two-dimensional semiconduc- 
tor with a 1.8 eV direct band gap,^°'^^ is an ideal material for vallej^ronics by 
valley-selective circular dichroism^, with ensuing valley polarization and valley 
Hall effect. 
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A non-equilibrium charge carrier imbalance between valleys is the key to 
creating valleytronic devices. The principal mechanism invoked here is circu- 
larly polarized optical excitation. ^'^ In this approach, the two valleys absorb 
left- and right-handed photons differently, a phenomenon referred to as circu- 
lar dichroism (CD). In order to obtain the valley-selective CD, it is essential 
then to break the inversion symmetry of the honeycomb lattice. In the case 
of graphene, it was suggested that by interacting graphene with a substrate 
such that the center of inversion can be obliterated, whereupon a gap opens up 
in each valley. ^^'^'^ This strategy, however, is quite challenging experimentally. 
And even in the eventual realization, there can only be a weak perturbation to 
graphene via the weak covalent coupling at large van der Waals separations. 

Since its first isolation, monolayer M0S2 has attracted immense attention. 
Many measurements have been performed to characterize the optical and 
transport properties of this material. In monolayer M0S2, two lay- 

ers of sulfur atoms in a two-dimensional hexagonal lattice are stacked over 
each other in an eclipsed fashion. Each Mo sits in the center of a trigonal pris- 
matic cage formed by 6 sulfur atoms (Fig. la). Quite remarkably in the context 
of current discussion, the natural stable structure of free-standing monolayer 
M0S2 is a honeycomb lattice with inequivalent bipartite coloring, breaking the 
inversion symmetry (see Fig. lb). Bulk M0S2 has an indirect band gap. In- 
terestingly, when thinned to the monolayer limit the material acquires direct 
band gaps located exactly at the corners of the Brillouin zone. Indeed, with its 
1.8 eV direct band gap and unique two-dimensional structure embracing the 
high-symmetry valleys, monolayer M0S2 has much to offer in the exploration 
of novel electronic and optielectronic devices and the associated physics. It is 
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quite natural to question whether it is possible to achieve valley-selective CD 
in this semiconducting atomic membrane, which will endow electrons in this 
material the valley degree of freedom, in addition to charge and spin that have 
been routinely explored in conventional device physics. 




Figure 1: The crystal structure of monolayer MoS2. (a) Coordination environment 
of Mo in the structure, (b) A top view of the lattice, emphasizing the connection 
to a honeycomb lattice. In our calculations, we use an optimized structure at the 
level of local density approximation in density functional theory (see text). The 
shaded region bounded by dashed lines corresponds to one primitive cell. The unit 
cell parameter is a = 3.12 A, and the vertical separation between sulfur layers is 
3.11 A. 



The key quantity to assess is, therefore, the k-revolved degree of optical 
polarization, ?7(k), between the top of the valence bands and the bottom of 
the conduction bands ^ 



r]{k,ujc 



|P^^(k)|2 + |P™(k)|2' ' ' 

This quantity is the difference between the absorption of left- and right-handed 
lights (±), normalized by total absorption, at each k-point and evaluated be- 
tween the top of the valence bands (f ) and the bottom of conduction bands (c). 



Note that the dependence on the transition energy, hucvO^) = ec(k) — e«(k), 
is imphcit through k. Here, the transition matrix element of circular polar- 
ization is V^{k.) = ;^[P™(k) ± iPy^{k.)]. The interband matrix elements, 
P'^''(k) = (ipcklpli^vk), are evaluated using the density functional perturbation 
theory (DFPT), within the local density approximation (LDA),^'' as imple- 
mented in VASP.^® Spin-orbit coupling is not included in our calculations. 
Briefly, a planewave basis set is employed at a cut-off energy 600 eV, and a 
total of 80 bands are included to ensure convergence of all computed quan- 
tities. A very dense k-point mesh (123 x 123 grid points) over the reducible 
two-dimensional full Brillouin zone is sampled in our calculations. 

As shown in Fig. 2a chiral absorption selectivity is indeed exact at K-t with 
1] = ±1. The contrast in the chiral absorptivities between the valleys owes 
its origin in the symmetry of both lattice and local atomic orbitals, which 
is quite different from the gapped graphene case^. A band state relevant to 
the optical excitation originates from local atomic states (or more generally, 
Wannier functions) bearing different orbital magnetic quantum numbers, /. In 
this case, the states at the top of the valence bands involve only d^2_y2 and 
dxy on Mo, and and Py states on S. At K+, the d-states on Mo hybridize 
as -^{dx2^y2 + idxy) {I = +2) to interact with -^{Px + iPy) (^ = +!)• At K_, 
the (i-states on Mo hybridize as ■^{dx2_y2 — id^^y) {I = —2) to interact with 
■^{Px~iPy) (^ = ~1)- At the bottom of the conduction band, only dz2 state on 
Mo {I = 0) is involved. These atomic orbitals form such linear combinations in 
accordance with the D^h point group symmetry and the lattice translational 
symmetry. 
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Figure 2: Valley-selective circular dichroism of monolayer MoS2. (a) Top valence 
band (blue) and bottom conduction band (pink). The center hexagon is the Brillouin 
zone color-coded by the degree of circular polarization, r](k), as defined in the text. 
The vector connecting K+ and K_ is perpendicular to Mo-S bond in the crystal 
structure in (a), (b) Schematic of phase winding on the M0S2 lattice that gives rise 
to the chiral optical selectivity. Left panel: the contribution to phase winding from 
the Bloch lattice phase, where r = ±1 is the valley index, and s = 1, 2 corresponding 
to the Mo and S sites (isospin) . Right panel: the phase winding under a three-fold 
rotation. The green axes indicate the rotation of local atomic coordinates that leads 
to the azimuth dissynchronization. 



The optical selection rule is rooted in the phase winding of the Bloch 
states under rotational symmetry, 3-fold rotation (C3) in this case. Given 
the symmetry-adapted linear combination orbitals, the azimuthal phase asso- 
ciated with the 3-fold rotation at K± is readily calculable, 
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'KK±)), 



(2) 



C'3|c(K±))=e^^2./3|^(K±)) 



(3) 



where v and c correspond to the valence and conduction band extrema, re- 
spectively. Notice, however, the phase winding associated with the rotation 
has two distinct contributions. The first comes from the Bloch phase shift in 
stepping from one lattice site to the next, as in the case of gapped graphene.'* 
The second phase factor arises as a consequence of dissynchronization of the 
azimuthal phase (associated with the magnetic quantum number, /) concom- 
mitant with rotation of the local atomic coordinates. These are schematically 
illustrated in Fig. 2b. 

Now the chiral optical selecvitity of the valleys can be deduced. The bottom 
of the conduction bands at the valleys, dominated by the / = d-states on 
Mo, bears an overall azimuthal quantum number ■m± = ±1, at K.±. At the 
top of the valence bands, m± = 0. Then for an optical transition at K-t, the 
angular momentum selection rule indicates that Am± = ±1, corresponding 
to the absorption of left- and right-handed photons. Therefore, our DFPT 
results concerning the close neighborhood of the valleys, which are the most 
important to the proposed optical valley polarization, are in fact ensured by 
the symmetry of the material. 

There is an important distinction, compared to gapped graphene, in the 
microscopic origin of chiral optical selection rule; that is, the selectivity arises 
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directly from the local relative azimuthal phase of the atomic orbitals, in con- 
trast to the sublattice-dependent Bloch phase winding in the case of gapped 
graphene/ The bonding in the electronic states in M0S2 across the gap is also 
considerably more complex, exhibiting richer possibilities of variation owing 
to the symmetry of local atomic states. Remarkably, the selectivity is nearly 
perfect over the entire valleys, and only changes sign rapidly across valley 
boundaries (Fig. 2a). This is to say, the entire valley K_|_ absorbs almost purely 
left-handed photons, whereas the entire valley K_ purely right-handed. The 
perfect valley-constrasting CD is very much conducive to optical polarization 
of the valleys. 

Now that we have established the valley-selective CD in monolayer M0S2, 
it is also interesting to look at the Berry curvature, f2„(k), which, if present, 
has crucial influence on the transport properties. Of note. Berry curvature 
enters into the semiclassical wavepacket dynamics via an anomalous veloc- 
ity perpendicular to the applied electric field (~ E x f2„(k)), in addition to 
the usual group velocity of Bloch bands. ^^'^^ The presence of non-vanishing 
Berry curvature is possible in the non-centrosymmetric honeycomb lattice.*^ In 
Fig.3, we plot the Berry curvature, fi„k,z = —'^I'^^idunk/ dkx\dun]i/ dky) along 
the K_ — r — path. Since the system has time-reversal symmetry and 
not inversion symmetry, fi„(k) is an odd function in k with generally non-zero 
values, as expected. The charge carriers' anomalous velocity acquires oppo- 
site signs in the two valleys, exactly canceling each other's contribution to 
the transverse current. At equilibrium, freestanding monolayer M0S2 will not 
exhibit anomalous Hall effect. Notice that at K-t, the Berry curvatures do not 
have the particle-hole symmetry. This is a clear indication that the physics of 
M0S2 cannot be fully captured by a minimalistic two-band model, as distinct 
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from the case of gapped graphene.^' 




Figure 3: Berry curvature, Qzi^), of bands across the band gap. The Berry cur- 
vature of only the states along the K„ — F — K+ path of the Brillouin zone are 
plotted. Notice that the value of Berry curvature is large for the conduction band 
at the zone center, where bands are degenerate. 



When valley polarization is induced by, say, valley-selective CD, only one 
valley has non-vanishing charge carrier population. This can then lead to the 
valley Hall effect.^ The Berry curvature across the band edges near K-t is 
most relevant to photo excited charge carriers. We see that at the band edges 
both conduction and valence bands display significant Berry curvature with 
opposite signs. Consequently, when electrons and holes are generated by a 
circularly polarized irradiation, both types of charge carriers have an intrinsic 
concurrent contribution to the Hall conductivity. Our results indicate that 
valley Hall effect in this material indeed warrants experimental assaying. In 
closing, we also would like to point out that M0S2 is but one of the many 
transition metal dichalcogenides, MX2, where M = Mo, W, and X = S, Se. 
They all have identical crystal structure and are similar electronically. These 

9 



compounds, once available in monolayer form, will provide a chemically rich 
family of semiconducting atomic membranes for exploring the valley physics, 
as we have demonstrated here for M0S2. 
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